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Abstract. Decoherence effect on quantum entanglement of two optical qubits in a 
lossy cavity interacting with a nonlinear medium (Kerr nonlinearity) is analyzed. The 
qubits are assumed to be initially in the maximally entangled states (Bell or Bell-like 
states) or the maximally entangled mixed states, on the example of Werner and Werner- 
like states. Two kinds of measures of the entanglement are considered: the concurrence 
to describe a decay of the entanglement of formation of the qubits, and the negativity to 
determine a decay of the entanglement cost under positive-partial-transpose-preserving 
operations. It is observed that the Kerr nonlinearity, in the discussed decoherence 
model, does not affect the entanglement of the qubits initially in the Bell or Werner 
states, although the evolution of the qubits can depend on this nonlinearity explicitly. 
However, it is shown that for the initial Bell-like state and the corresponding Werner- 
like state, the loss of the entanglement can be periodically reduced by inserting the 
Kerr nonlinearity in the lossy cavity. Moreover, the relativity of the entanglement 
measures is demonstrated, to our knowledge for the first time, as a result of a physical 
process. 
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1. Introduction 



Decoherence, resulting from the unavoidable and irreversible coupling of a quantum 
system to its environment, turns a correlated quantum state of the system into a classical 
statistical mixture p. Decoherence, causing usually a loss of quantum entanglement, is 
one of the major limitations of practical capabilities of quantum computers [2]. Thus, 
the analysis of the dynamics of entangled quantum two-level systems (qubits) coupled 
to the environment, represented by a thermal reservoir, is of a particular importance. 
In this paper, we will study the loss of the entanglement due to a dissipative nonlinear 
interaction of two optical qubits, which are implemented by superpositions of vacuum 
and single-photon states of two cavity modes, and assumed to be initially in the 
maximally entangled states (MESs) or the maximally entangled mixed states (MEMSs). 

It is a well accepted fact that there is no unique way to quantify mixed-state 
entanglement and thus various measures with different operational interpretations have 
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been proposed to describe different aspects of the entanglement. We will apply the 
concurrence [3], a measure related to the entanglement of formation and the 
negativity 0101111, a measure corresponding to an operation-limited entanglement cost 
[Slinj. We have chosen these particular measures as they are similar to each other from a 
physical point of view and, moreover, can easily be calculated contrary to other measures 
including the entanglement of distillation or the relative entropy of entanglement. 

Our analysis is related to a new regime of quantum nonlinear optics involving 
highly-efficient nonlinear interactions between very weak optical fields, which has 
been recently demonstrated experimentally in, e.g., dense atomic media by using an 
electromagnetically induced transparency (EIT) to resonantly enhance nonlinearities 
(for a review see JHl)- In particular, observation of giant Kerr nonlinearities has been 
predicted [TT] and first measured in an ultracold gas of sodium atoms to be ~ 10^ greater 
than those in the conventional optical materials [12J. Physical realizations of a Kerr 
nonlinear cavity enabling strong interaction of photons was suggested by Imamoglu et 
al. [m Unj and then studied by others [TU [13 HEl HZl HE] • Motivated by these advances, 
there is an increasing interest to apply the Kerr nonlinearities for quantum information 
purposes (TH], including the problem of generation of highly-entangled states (see [THtl^ 
and references therein). Nevertheless, the effects of decoherence on the entanglement of 
fields interacting via the Kerr nonlinearity have not been discussed in greater detail yet. 

The paper is organized as follows. In section 2, we define the entanglement measures 
to be used in our description of decoherence. The model and its solution for two optical 
qubits in a lossy nonlinear cavity are presented in section 3. The main results concerning 
the decoherence of the qubits being initially in the maximally entangled states and the 
maximally entangled mixed states are presented in sections 4 and 5, respectively. A 
physical implementation of the model and discussion of the results are given in section 6. 

2. Entanglement measures 

We will apply two measures of entanglement to analyze the effect of decoherence on the 
entangled qubit states. The first measure is the concurrence defined for two qubits as |3] 

4 

C(p) = max{2 max Aj — Aj, 0} (1) 

* i=i 

where Aj are the square roots of the eigenvalues of the matrix p(o"iy ® d'2y)p*{d-iy Cg) a2y) 
with ajy being the Pauli spin matrix of the jth qubit and the asterisk denotes complex 
conjugation. The entanglement of formation, Ep{p), which characterizes the amount of 
entanglement necessary to create the entangled state |3j, is for two qubits given by a 
simple monotonic function of the concurrence jH] 

EAP)=H{1[1 + ^1-Cim (2) 

where H{x} = — xlogga; — (1 — a;) log2(l — x) is the binary entropy. The entanglement 
of formation and, equivalently, the concurrence vanish for an unentangled state, and 
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are equal to one for a maximally entangled state. A measure associated with the 
entanglement of formation is the entanglement cost defined as |5 lining oo-£'f(p'^"')/'^ 
which, in general, is quite difficult to calculate. Thus, for simplicity, we will describe the 
entanglement cost limited to a special class of operations to be specified in the following. 

Another useful measure of the entanglement is the negativity [HI IZj, which 
corresponds to a quantitative version of the Peres- Horodecki criterion We adopt 
here the following definition 

Af{p) = 2 max(0, - ^ fij) (3) 
j 

where the sum is taken over the negative eigenvalues fij of the partial transpose p^-* 
of the density matrix p of the system. For two-qubit pure or mixed states, the sum 
in (jni) can be skipped as p^-* has at most one negative eigenvalue [21]. The negativity 
satisfies the standard conditions for a useful measure of the entanglement [221 123 ■ For 
two-qubit states, the negativity, defined by Q, becomes 1 for a MES and vanishes 
for an unentangled state, the same as the concurrence. Recently, Audenaert et al. 
jH] and supplementary Ishizaka P have provided an operational interpretation of the 
logarithmic negativity, defined by |2S1 

E^ip) = \og,[Nip) + ll (4) 

as a measure of the entanglement cost for the exact preparation of a two-qubit quantum 
state p under quantum operations preserving the positivity of the partial transpose 
(PPT). 

For an arbitrary two-qubit pure state 

I*) = CoolOO) + CoilOl) + CiollO) + cnlll) (5) 

with the normalized complex amplitudes Cij, we have the concurrence and negativity 
equal to each other and given by a simple formula 

Afis, = Cis, = 2|cooCii - CoiCiol. (6) 

However, for qubits in a mixed state, the entanglement measures are usually different. 
In general, the inequality A/'(/3) < C (p) holds for an arbitrary two-qubit state p as first 
observed by numerical investigation by Eisert and Plenio |2j and Zyczkowski |21], and 
then proved by Verstraete et al. j2SI, 

Eisert and Plenio [Zj raised an intriguing problem of the relativity of entanglement 
measures: if according to one measure of the entanglement the state pi is more entangled 
than p2 then does it imply that pi is also more entangled than p2 according to another 
entanglement measure? By Monte Carlo simulation of thousands of two-qubit states, it 
was observed that indeed the condition 

c{p,)<cip2)^^^iPl)<^^iP2) (7) 

can be violated by some states pi and p2 jH 1211 1201 • It should be stressed that this odd 
looking property is physically sound since such incomparable states pi and p2 cannot 
be transformed into each other with unit efficiency by any local quantum operations 



Decoherence of two maximally entangled qubits 



4 



and classical communication (LQCC). In general, all good asymptotic entanglement 
measures are either identical or put different orderings of quantum states as implied by 
the requirements of equivalence and continuity of the measures on pure states |2Zj. Thus, 
by comparing various entanglement measures defined to examine different methods of 
the entanglement preparation and/or its use, one indeed can arrive at the problem of 
different state orderings imposed by the measures. The only alternative to avoid the 
state-ordering ambiguity is to declare one entanglement measure for mixed states as the 
unique one, but this would preclude us from examining the problems of how to prepare 
the entanglement and how to make use of it ^|. Here, we will give simple analytical 
examples of states differently ordered by the concurrence and negativity, thus we will 
explicitly demonstrate the relativity of these entanglement measures. 



3. Model and its solution 



Decoherence effects on optical modes (qubits) in a lossy nonlinear cavity can be described 
by a model of coupled dissipative nonlinear oscillators represented by the following 
prototype Hamiltonian ]29^t 

H = Ho + H^^ + Hi (8) 

where 

N N 

Ho =/^E^.4% + ^EE4^'^(^V&?, (9) 

j=l k j=l 

N 

-^NL = ^ E Xijaldid]dj, (10) 

N 

Hi =^EE[^??«l^i'^ + h.c.] (11) 

and dj is the annihilation operator for the jth system oscillator at the frequency (jj; is 
the annihilation operator for the fcth oscillator in the jth reservoir at the frequency ^^P'-, 
Xij are the nonlinear self-coupling (for i = j) and cross- coupling (for i ^ j) constants 
proportional to a third-order susceptibility of the Kerr nonlinear medium (see Sect. 6), 
and g^^^ are the coupling constants of the reservoir oscillators. Dissipation of the system 
is modelled by its coupling to reservoirs of oscillators as described by Hamiltonian Hi. 
The evolution of the dissipative system under the Markov approximation is governed 
by the following master equation for the reduced density operator p in the interaction 
picture 

—p = TT-[^NL, p] -h E -^{'^ji'^^jP^j - ^j^P - pO'ja'^^ 



+ {uj + l)(2ajpd] — a]ajp — pa]aj)} (12) 

X Hamiltonian J?nl is sometimes defined in the normal-ordered form of and aj. However, such 
Hamiltonian differs from ours only in terms proportional to Xjj'^\o,j^ which can be incorporated in the 
free Hamiltonian /fg, so this modification does not effect the entanglement. 
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where nj are the mean thermal occupation numbers and 7j are the damping constants, 
which will be assumed the same, 71 = 72 = 7, in the next sections. With the help of 
a disentanghng theorem for SU(1,1) in thermofield dynamics notation, Chaturvedi and 
Srinivasan |^ have found a general solution of the master equation (fT^ both for the 
quiet [uj = 0) and noisy {nj > 0) reservoirs. By confining our analysis to the case of 
only two oscillators (A^ = 2) coupled to the quiet reservoirs, the Chaturvedi-Srinivasan 
solution for the density matrix elements (mi, m2|p(t)|ni, 77-2) in the photon-number basis 
can be written as 



{mi,m2\pit)\ni,n2) = J2 H ^i^2("^i + Pi, "^2 + P2|p(0)|ni + pi, ^2 + P2) (13) 

pi=0 P2=0 



where 



Rj = Rj {mj , Uj , pj ) 




- exp(— Xjt)] 




X exp |i(xji + Xj2){mj - nj)t - [xj {mj + Uj + 1) - 7j]^| (14) 

with Xj = 7j + 2i[xji{mi — ui) + Xj2(^2 — ""-2)], and (^^^ are binomial coefficients. In our 
scheme, qubits can be represented by the single-cavity modes restricted in the Hilbert 
space spanned by vacuum and single-photon states (see, e.g., ^Q). Then, for the qubit 
states, solution (fT^ simplifies to the summations over pi,p2 = 0, 1 only. 

By assuming no dissipation (71 = 72 = 0) in our system, the evolution is governed 
by the unitary operator exp{—iH^^i^t/h). It implies that, for the two qubits initially in 
a pure state the concurrence and negativity evolve periodically as follows: 

Ar*(7 = 0,t) = C.,(7 = 0,t) =2|exp(2zxi2t)coo(0)cii(0)-coi(0)cio(0)| (15) 

depending on the cross-coupling X12 but not on the self-coupling constants Xii and X22- 
One can observe that the evolution of the qubits in the nonlinear medium can lead to 
a periodical generation of entangled states. Even for the initial separable state 

1^(0)) = (t/i|0)i + d2\l)l) ® (rf3|0)2 + d^\l)2), (16) 

where + |(i2p = + = 1 and none of the amplitudes di is zero, the 
concurrence and negativity periodically become positive 

7Vvt(7 = 0,t) = Cvi,(7 = 0,t) = 4\did2d3d4sm{xi2t)\ (17) 

which corresponds to the entanglement of up to -f^{|[l + y^l — 16|(ii(i2(i3(i4P]} ebits. In 
particular, the initial state (fTBj) with all the amplitudes equal to 1/v^, i.e., 

mo,)^Mi^«M,.|,,,>, ,,,, 

evolves into a maximally entangled state, defined below by ()26|) . having exactly 1 ebit 
at the evolution moments t = (1 + 2n)TT / {2x12) {n = 0, 1, ...). Nevertheless, the MESs 
are not generated if our system is subjected to dissipation. 
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4. Decoherence of the maximally entangled states 

Let us assume that two qubits are initially in the Bell states 

l^±) = ^(|01)±|10)) (19) 

which evolve in the dissipative nonlinear cavity into the following mixed states 

p^±(t) = ^{2(1 - g)mm + ^?(|01)(01| + |10)(10|) 

± ^(e*(^i-^2)t|oi) (io| + h.c.)} (20) 

where g = e"'''* and Xi = Xu- The evolution is independent of the nonlinear cross- 
coupling Xi2 but depends on the self-couplings xi and X2- We find that the concurrence 
is simply given by 

CS)=9 (21) 
and the negativity is 



±) = -=(|00)±|ll)) (23) 



Ar4t) = ^2g^-2g + l+g-l (22) 

being independent of the sign in ()19p. As implied by the form of the density matrices 
(I2()|l . the entanglement measures are independent of any Kerr couplings. On the other 
hand, the initial Bell states 

1 

V2 

evolve in our lossy system into 

p<^±(t) = ^{(2 ~2g + /)|00)(00| + (1 - ^7)^7(|01)(01| + |10)(10|) 

± ^(e^(^^+2»2+^2)*|00)(ll| + h.c.) + /|11)(11|}. (24) 

Contrary to p^±(t), the density matrices p^±{t) depend on the cross-coupling between 
the qubits. We find that the concurrence and negativity are the same for any evolution 
times and any sign in ()23|) as given by 

C^{t) = M^{t) = g^ (25) 

in contrast to C^{t) and J\f^(t), given by pij) and (j221), respectively, which are the same 
at t = and t = oo only. There are the following important properties of the discussed 
entanglement decays. First, the concurrence (negativity) for the initial Bell states \ip±) 
decays slower (faster) than that for \4>±), as it holds C^(t) > C^{t) and Af^(t) < N^{t) 
for any damping constants 7^ > (fc = 1, 2) and any moments of time < t < 00. Thus, 
we provide an explicit example of states violating condition ((Tj). Second, contrary to the 
density matrices, the entanglement measures are independent of the nonlinear couplings 
for the initial Bell states ()19|) and (j2SI)- Thus, decoherence-free evolution in the nonlinear 
cavity does not change the entanglement, i.e., 6*^(7 = 0,t) = ^^^(7 = 0,t) = 1. Now 
we will give an example of a maximally entangled two-qubit state evolving in the Kerr 
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(27) 



medium in such a way that the entanglement depends on the cross-couphng Xi2- Let us 
analyze the following initial state 

\^) = i(|00) + |01) + |10) - 111)) = i=(|0, +) + |1, -)) (26) 

where |±) = (|0) ± |l))/v^. The state \^) is a MES since its concurrence and negativity 
are equal to one. For brevity, we neglect the self-couplings, Xi = X2 = 0, which do not 
affect the entanglement. Then the initial pure state \if) evolves in the Kerr medium into 
the mixed state described by 

(2 -(7)2 -fg 
h*^ g{2-g) g -fg^'^ 
h*V9 9 9(^-9) -f9'^' 

-1*9 -r9"' -1*9'/' 9' 

given, as usual, in the computational basis {|00), |01), |10), |11)}, and h = {•yfg — 
2^X12)7(7 — 2^X12), g = exp(— 7t), and / = exp(2ixi2t). If we moreover assume no 
losses in the nonlinear cavity (7 = 0), then the evolution of the initial state \ip) results 
in the entanglement oscillations described simply by 

C^(7 = 0,t) =Ar^(7 = 0,t) = |cos(xi2t)|, (28) 

which is in contrast to the case of the initial Bell states and \ip±), which evolve 
without changing their entanglements. The aperiodic decay of the entanglement for the 
density matrix (j27jl occurs only if there is no interaction between the qubits and then 
it is described by 

1 

2' 



C^iXi2 = 0,t) = -gil + g), (29) 



KiXi2 = 0, t) = - 4x + 1 + ^ - 1 (30) 

where x = g{l — g)/2. For nonzero damping and cross- coupling parameters, both the 
concurrence C<^(t) and the negativity N^{t) exhibit decaying oscillations, as shown by 
curves (c) in figures 1 and 2, respectively. The expressions for C^(t) and A/'^(t) are quite 
lengthy in general case of non-zero X12 and 7, so we do not present them here. Instead 
we give approximate formulas for the envelope function of the concurrence 



C^rW ^ |[\/^ - + ^21^ -2y){x + y-z))+g-l] (31) 



where x = 27 - 14^ + 3^^ y = v/T59^T29^T37^^^^3^, and z = ^{x + yf - 9y^, 
and for the envelope of the negativity 

Ar;-(t) ^ 1 (2Re{/v + t3il-g)gV3^ - (2 - g)' - g^ (32) 

where v = 8g^ - 18g^ - 93^^ + 324c/3 - 273^^ ^ igQg _ 54 ^nd w = 116/ - 316^^ + 
297 g"^ + 930g^ — blbg"^ + Q2Ag + 16. Equation (jH^ was derived by assuming only that the 
cross-coupling Xi2 is much stronger than the damping constant 7, which implies that the 
function h in the density matrix (jTfjl approaches unity. The envelope functions, given 
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by (jnH) and ((221), depicted by curves (e) in figures 1 and 2, respectively. Another 
simpler but far less accurate approximation of the negativity envelope function can be 
given by 



which was obtained by using the general properties of the eigenvalues /Xj of the partial 
transpose ff^it) of density matrix (j2ZI), including Z^i/^j = 1 and Wi^i = det p^^{t), and 
observing that there exist two eigenvalues /ij summing up approximately to zero. It is 
worth noting that the envelope functions, the same as those given by and (jSS), are 
for the system initially in the separable state given by ()18p. for which p{t) has the form of 
(07|) but with the functions h and / modified as follows h = ['y{2+fg) — 2ixi2]/ (7— 2ixi2) 
and / = — exp(2ixi2t). Note also that the envelope functions are independent 

of the cross-coupling X12 under assumption X12 ^ 7 but, even in this regime, the 
period of entanglement oscillations is a function of Xi2- A closer comparison of the 
entanglement for the Kerr interacting and non-interacting qubits in the lossy cavity 
leads us to the following inequalities C™^(t„) = C,^(xi2 > 0,t„) > C^{xi2 = 0,t„) and 
A/'™^(t„) = A/'^(xi2 > 0,tn) > A/'^(xi2 = 0,t„) valid for the moments of time equal 
to tn = nrc /xi2 for n = 1, ■ ■ ■• By comparing the entanglement measures for the all 
analyzed MESs (see figures 1 and 2) we can finally conclude that 



where the equalities hold for the nonzero damping constant 7 at the evolution moments 
t = and t = 00, while for 7 = at any times t. Inequalities except those for 

C^''{Xi2 > 0,t) and A/'^'^^(xi2 <|t7,t), can be proved analytically, while the remaining 
inequalities were checked numerically for a large class of parameters. Note that for 
small values of X12 in comparison to 7 it holds J\fip{xi2,t) < A/'^(t), nevertheless the 
last inequality in (jH^j) is satisfied even if X12 ~ 7, and more pronounced for X12 ^ 7 
(see figure 2), which is the condition assumed in the derivation of (p?T|) - (p?H|) . Obviously, 
inequalities corresponding to (jH^ hold for the entanglement of formation, Ep{t), and 
those corresponding to are also valid for the PPT-entanglement cost, EN{t). The 
main conclusion is the following physical interpretation of inequalities (jMjl - ljH^jl : by 
enabling Kerr interactions between the qubits initially in the Bell-like state, given by 
PHjl . the loss of the entanglement can be periodically reduced. 

5. Decoherence of the maximally entangled mixed states 

We will analyze the decoherence process of the initially maximally entangled mixed 
states of two qubits [HS ESI IHH] on the example of the Werner states [SEj defined to be 




(33) 



C^t) > C-^(xi2 > 0,t) > C^{xi2 = 0,t)> C^{t), 
M^it) < M^{xi2 = 0,t) < M^it) < Ar;-(xi2 4^1. t) 



(34) 
(35) 



(l/3<p<l): 



Pp^±(0) =p|^±)(V'±l + 




h ® /2, 



(36) 
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Figure 1. Decay of the concurrence for the initial Bell and BcU-Uke states: (a) 
C^{t), (b) C^{t), (c) C^{x'i2,t), (d) C^{xi2 = 0,i) (dashed curve), and (e) ^^-(x'la.i) 
(dotted curve) for cross-couphng constant Xi2 = 20 rad MHz and damping constant 
7 = 4 rad MHz. 




Figure 2. Decay of the negativity for the same Bell(-Uke) states and interactions as in 
figure 1: (a) AA^(t), (b) AA^(i), (c) K{x'i2,t), (d) NUxi2 = 0,t), and (e) AA°-(xi2,i)- 
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p=0.9 p=3/4 




time xio-°s time xiq-'s 



Figure 3. Decay of the concurrence for the initial Werner(-hke) states: (a) Cp^(i), 
(b) Cp^{t), (c) Cp^{xi2^t), (d) Cp^{xi2 = 0,t), and (e) C™''(xi2,0 for various values 
of parameter p; 7 and X12 the same as in figure 1. 

Pp0±(O) =p|0±)(0±| + ®/2 (37) 

where \ip±) and |0±) are given by (fT^ and (|^. respectively, and Ji,2 are the identity 2x2 
matrices. Thus, the Werner states are mixtures of a MES (Bell state) and the maximally 
mixed state, given by Ji ® I2, which can be interpreted as an equal incoherent mixture 
of the four Bell states. It is worth mentioning that the standard two-qubit Werner 
state is defined as pp^_(0) only This state, given in terms of the singlet state 

lip-), is invariant if both qubits are subjected to same unitary transformation, U ®U . 
Nevertheless, by ignoring the U ® U invariance but keeping the same entanglement 
properties, the standard Werner state is often generalized (see, e.g., 123 IHll Ell) to 
include mixtures of any MESs, as given by (jHE)) and (jSZ|)- Following this convention, we 
will apply the generalized definitions of Werner states in our study. 

It is easy to show that the concurrences and negativities of the Werner states are 
the same and given by 

Cp^(O) = C^p0(O) = KM = KM = (3P - l)/2. (38) 
The Werner states can be considered the MEMSs since their degree of entanglement 
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Figure 4. Decay of the negativity for the same Werner(-Hke) states and interactions 
as in figure 3: (a) Afp^{t), (b) Afp^{t), (c) Up^{x'i2^t), (d) ^fpv{Xi2 = 0,t), and (e) 

Aer(xi2,t)- 



cannot be increased by any unitary operations j32] and they have the maximum of 
entanglement for a given hnear entropy (and vice versa) jHSl- In a special case of p = 1, 
the Werner states go over into the MESs. The evolution of Pp^{t) for the initial Werner 
state (jHBjl in the lossy nonlinear cavity is described by 

Pp^t) = \m - 9? - /p]|00)(00| + g\l - 

± 2c/p[e^(^i-^2)t|oi)^io| + e-^^^-'^^)* 1 10) (01 1] 

+ (7[2-(7(l-p)](|01)(01| + |10)(10|)}, (39) 

being independent of the cross-coupling Xi2-, which implies that a monotonical decrease 
of the entanglement occurs according to 



C,^{t) = max{0, gp - g\ (1 - g){l - p) + ^'^^ }, 



4 



M,.^{t) = max{0, ^(l-gy + g^p^ - 9!0_Pl _ (l _ ^)}. (40) 
In a special case of p = 1, the above formulas simplify to (PT|) and ((221), respectively. 
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On the other hand, the evolution of Pp^{t) from the initial Werner state (j37|) reads as 

Pp<t>±{t) = ^{(2 -2g + xp)|00)(00| + {g - Xp)(|01)(01| + |10)(10|) 

±p(/|00)(ll| + r|ll)(00|) + x,|ll)(ll|} (41) 

where Xp = {1 + p)g'^ /2 and / = 5'exp[i(xi + 2xi2 + X2)i]- Hence, the time evolution 
explicitly depends on the cross-coupling Xi2i but in such a way that the concurrence 
and negativity exhibit the same monotonous decrease independent of Xi2 as follows: 

Cp^{t) = Afp^it) = max{0, ^[g{l+p) - 2{l-p)]} (42) 

In a special case of p = 1, equation reduces to g^ in agreement with pKjl . Note 
that the subscript ± in Cp^, Np^, Cp^, and Np^ has been omitted as the functions are 
independent of the sign in (j36|) and pTj) . 

As the last example, let us assume qubits to be initially in the Werner-like state 
defined by (1/3 < p < 1): 

Pp^(0)=p|<^)(<^| + ^/i®/2 (43) 

in terms the MES given by ()26p. The concurrence and negativity for ()43p are equal to 
Cp(p(0) = A/^(^(0) = (3p — l)/2 being the same as for the other Werner states. However, 
its evolution essentially differs from Pp^±(t) and Pp^±(t) by exhibiting oscillations of the 
entanglement. In detail, it is described by the density matrix elements 

= P^'~~'''\Mt)h (44) 
given in terms of (j2Z|), but with the off-diagonal terms multiplied by p as 6ij stands for 
Kronecker delta. In a special case of the lossless nonlinear cavity, the entanglement of 
the state Pp^(7 = 0,t) evolves periodically as follows: 

Cp^(7 = 0,t)=Afp^{-f = 0,t) = imax{0,p(2|cos(xi2t)| + 1) - 1} (45) 

which is opposite to the time-independent evolution of the other Werner states, viz. 
PpV±(7 = 0)^) = Pp<f>±{l = 0,t) = const. One can conclude from (jl^ (see also figures 3 
and 4) that by decreasing parameter p, the entanglement and the time intervals in which 
the states are entangled decrease. For the dissipative nonlinear cavity, the entanglement 
corresponding to the evolution of Ppipit) exhibits decaying oscillations shown by curves 
(c) in figures 3 and 4. As in the former section, we are mainly interested in the envelopes 
of these oscillations. In the special case of p = 1, when the initial Werner- like state goes 
over into the Bell-like state, the concurrence and negativity envelopes are given by 
and (1221), respectively. By assuming Xi2 3> 7, an approximate formula for the 
p-dependent envelopes of the concurrence can be given by 

CpT W - f {O' ^ (V^P + 4Pv^ - 2^Xp - 2p^) +g + p-2} (46) 

in terms of a;p = 3G^ + 2Gp + llp'^ and ?/p = 3G^ + {10 + 9p) + G {3 + Up) +p{9 + 16p) 
where G = 2 — g. Note that for p = 1 is another approximate formula of 

the concurrence envelope for the initial MES \ip), but leading to a slightly worse 
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Figure 5. Four-level atomic system for the resonantly enhanced Kerr nonlinearity in 
the Schmidt-Imamogiu scheme. 

approximation than that given by (j31|) . For brevity, the lengthy formula for the p 
-dependent negativity envelope, Mp^it), generalizing equation (jH^ . is not presented 
explicitly here although was used for plotting the envelope curves (c) in figure 4. By 
analyzing figures 3 and 4, we conclude that 

■AA;7(tn) = Mp^ixi2 > 0, t„) > ^fp^iXl2 = 0, t„) (47) 

at moments of time in ~ nir/xu = 1,2, ■■■), which means that the decay of 
entanglement of the initially Werner-like state in a lossy cavity, can be periodically 
retarded by inserting the Kerr nonlinearity in the cavity. 

6. Discussion and conclusions 

Let us finally address the question whether the interactions studied in this paper can 
be experimentally observable. As mentioned in Introduction, the conditions assumed 
in the paper of the strong Kerr interaction at low light intensities can be satisfied, e.g., 
for the EIT schemes as studied theoretically [TH [T31 Cll UHl HEl HI] and confirmed 
experimentally E]- Schmidt and Imamoglu ^1] have proposed a renown EIT 
scheme where a low density cloud of cold atoms with the four-level structure, shown 
in figure 5, exhibits giant resonantly-enhanced nonlinear cross-coupling with vanishing 
linear susceptibilities at low intensities. In the scheme, atoms are placed in a cavity (or 
double cavity) tuned to two frequencies: ui of the mode ai resonant with the transition 
|1) |3), and U2 of the mode 02 detuned by Auj2 of the transition |2) ^ |4). The 
EIT effect for the modes ai and 02 is induced by a classical coupling field of frequency 
Uc resonant with the transition |2) ^ |3). By assuming that \u!i — uj2\ ^ A.U2, no 
nonlinear self-coupling will occur in the system. Nevertheless, as shown in the preceding 
sections, only the cross-coupling changes the entanglement evolution. By adiabatically 
eliminating all the atomic levels, Schmidt and Imamoglu have found the real part 
of the resonantly-enhanced third-order nonlinear susceptibility, Re(x*^^-'), to be given 
by 1/^13 P 1/^24 P'^at(2eo^^f2c^co'2Kav)~^, where fXij is the electric dipole matrix element 
between the states \i) and \ j), riat is the total number of atoms contained in the cavity 
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of volume Kav, is the coupling- field Rabi frequency, and eo is the permittivity of 
free space. With the help of expression for Re(x^^''), it is easy to show that the Kerr 
nonlinear cross-coupling is given by fSl 

2^^^ ^ Q^Au, ^^^^ 



where gij = fiij ^ Ui / {2ehVcstv) is the coupling coefficient between the atoms and the 
cavity mode Oj of frequency uji. It is worth stressing that the above formulas for Re(x*-^^) 
and Xi2 are valid under condition that [(yfispnat/f^c < ^ required by the applied adiabatic 
elimination procedure jHZI. The EIT enhanced Kerr-coupling constants for the Schmidt- 
Imamoglu scheme can be estimated moderately as ~ 0.2 rad MHz 15^ or, by putting 
the stringent limit on the required cavity parameters |HZ], as ~ 100 rad MHz In 
our numerical analysis, we have chosen Xi2 = 20 rad MHz. A typical cavity decay rate 
obtainable in current experiments is of the order ~ 4 rad MHz, which is 5 times smaller 
than the value of Xi2 chosen for plotting figures 1-4. This estimation is less stringent 
than that given in Ref. |l3j. It is worth noting that the EIT leads to remarkable light- 
speed reduction ^l2], which enables reduction of the cavity decay rate in the Schmidt- 
Imamoglu setup with the same finesse mirrors. We have analyzed decays within times 
< 80ns, which are fairly shorter than dephasing time (~ 9/is) for atom cloud measured 
in the Hau et al. experiment ^2] but longer, for the obvious reasons, than the evolution 
times (~ 8ns) in the quantum non-demolition scheme of Duan et al. [T3] . 

In conclusion, we have analyzed evolution of two optical modes in qubit states 
interacting via a Kerr nonlinearity in a lossy cavity modelled by dissipative coupled 
nonlinear oscillators being initially in the maximally entangled pure or mixed states. 
We have found that for the initial Bell {\ip±) and \<f)±)) or Bell-like {\^p)) states, the 
decay of the concurrence, or equivalently the entanglement of formation, is the slowest 
for l-^-t) and the fastest for |0±), while the decay of the negativity, or equivalently 
the PPT-entanglement cost, is the slowest for \(p) (if the nonlinearity parameter is 
much greater than the damping constants) and the fastest for Thus, we have 

provided simple analytical examples of states differently ordered by the concurrence and 
negativity. These seemingly inconsistent results are physically meaningful as discussed 
in, e.g., Refs. jH EH EE! and proved in general terms by Virmani and Plenio |27j . 
Nevertheless, to our knowledge, our analysis is the first demonstration of the relativity 
of the entanglement measures as a result of a physical process. Moreover, we have also 
studied decoherence of the initial maximally entangled mixed states on the example of 
three kinds of Werner(-like) states as related to the different Bell(-like) states \ip±), \4>±) 
and \ip). Our analytical and numerical results show the differences and similarities of 
the negativity and concurrence decays of the Werner(-like) states in comparison to the 
Bell(-like) states. 

We have demonstrated that by inserting medium with the Kerr nonlinearity, 
described by Hamiltonian (|TO|l . into the lossy cavity, evolution of the initial Bell states 
\ip±) or and the corresponding Werner states is changed but in such a way that the 
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entanglement decays in the same manner as without the nonhnear medium. However, 
if the qubits placed in a lossy cavity are initially in the BcU-like state \ip) or the 
corresponding Werner-like state, the loss of entanglement can be periodically delayed 
(partially recovered) by inserting medium with the Kerr nonlinearity. 
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